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The magnon thermal conductivity /-c mag of the spin ladders in Sri4Cu24- x Zn x 04i has been investi- 
gated at low doping levels x = 0, 0.125, 0.25, 0.5 and 0.75. The Zn-impurities generate nonmagnetic 
defects which define an upper limit for Z mag and therefore allow a clear-cut relation between Z mag 
and K mag to be established independently of any model. Over a large temperature range we observe 
a progressive suppression of K mag with increasing Zn-content and find in particular that with respect 
to pure Sri4Cu2404i « mag is strongly suppressed even in the case of tiny impurity densities where 

In 



< 374 A. This shows unambiguously that large In 



3000 A which have been reported for 



Sri4Cu2404i and LasCagCr^CUi on basis of a kinetic model are in the correct order of magnitude. 



I. INTRODUCTION 

The transport of energy and magnetization 
in low dimensional quantum magnets is sub- 
ject to numerous experimental and theoretical 
in-.-cctigqtijn; , 1 - 2 - 3 - 4 - 5 - 6 - 7 - 8 - 9 - 10 - 11 - 12 - 13 - 14 - 15 - 16 - 17 - 18 - 19 - 20 - 21 - 22 
with the nature of scattering and dissipation of thermal 
and magnetic currents at the center of interest. In 
particular, quasi one-dimensional (ID) systems are 
being intensely studied since properties starkly different 
from classical systems are expected: ballistic magnetic 
heat transport, i.e. a finite thermal Drude weight Dth, is 
predicted in integrable models like the XXZ Heisenberg- 
chain with S = i/2, 10 , 12 i 13 , 14 i 15 , 24 Experimentally, for 
the 5=1/2 two-leg spin ladder systems Sri4Cu2404i 
and La5CagCu2404i exceptionally large values of the 
magnetic thermal conductivity K mag are found (up to 
~ 100 Wm _1 K _1 even at room temperature) i2*ii^ In 
principle, these data should provide the unique possi- 
bility of extracting the magnetic mean free path ^ mag , 
hence obtaining information about magnetic scattering 
processes in such low-dimensional spin systems for the 
first time. It turned out, however, that the extracted 
values of ^ mag are strongly model dependent: on one 
hand, the analysis within the framework of a simple ki- 
netic model yielded values up to i mag » 3000 Ai£ For the 
complicated and supposedly rather imperfect materials 
studied such values are surprisingly large and might arise 
from oversimplifications inherent to the model. On the 
other hand, much smaller values were obtained in a more 
sophisticated analysis which is based on a finite D t h~ 
However, the finiteness of Z? t h in these systems is subject 
of a highly controversial debate since both ballistic 
and non-ballistic heat transport has been predicted for 
non-integrable S = 1/2 models such as these two- leg 
spin-ladders or frustrated chainsi 11 ' 14 ! 15 ! 16 ' 19 ! 20 ' 21 ! 23 . 



material and thereby impose an upper limit for i mag in 
a well-defined manner. We achieved this by substitut- 
ing small amounts of nonmagnetic Zn 2+ for the mag- 
netic Cu 2+ (with S = 1/2). K mag of Sri4Cu24- x Zn x 04i 
(a; < 0.75) is progressively suppressed with increasing Zn- 
^ojjtgpt and, in particular, already at tiny impurity densi- 
ties yielding Z mag < 374 A this suppression is strong with 
respect to the pure material (exceeding a factor of ~ 5). 
Employing the simple kinetic model which had been used 
in previous studies^ we even find Z mag « dzn-Zn at low 
temperatures, where dzn-Zn is the nominal mean dis- 
tance between Zn-ions along a ladder. This unambigu- 
ously shows that the large Z mag « 3000 A which have 
been reported for Sri4Cu2404i and La5CagCu2404i are 
in the correct order of magnitude. 



II. EXPERIMENTAL 

We have grown single crystals of Sri4Cu24- x Zn x 04i 
(a; = 0, 0.125, 0.25, 0.5, 0.75) by the travelling sol- 
vent floating zone method^ 7 - The thermal conductivity 
as a function of temperature T of these crystals has been 
measured with a focus on the crystallographic c-direction 
(k c ), which is parallel to the orientation of the ladders^ 
The measurements in the range 7-300 K have been per- 
formed with a standard four probe techniquei 2 * 2 ^ 9 . 



III. RESULTS AND ANALYSIS 

Fig. n a) shows the temperature dependence of the 
thermal conductivity of the undoped material, i.e. 
Sri4Cu24- x Zn x 04i at x = along all three crystal axes. 
In this material the thermal conductivity along the lad- 
der direction, n c (T), frequently displays a double peak 



In this article we present an investigation of K mag and structure 
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The low-temperature peak (~ 20 K) is 



^ mag of Sri 4 Cu 2 404i 



which is independent of any model. 
The basic idea is to generate nonmagnetic defects in the 



observed along all crystal directions and arises due to the 
usual phonon heat conduction. The peak at higher tern- 
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FIG. 1: Thermal conductivity of Sri4Cu24-xZn x 04i (x = 
0, 0.125, 0.25, 0.5, 0.75) as a function of temperature, a) 
Anisotropic thermal conducti vity along a-, b-, and c-axes (n a , 
Kb, k c ) at x = (cf. Ref. 1 21-51) . b) k c at x > 0. Solid 
lines represent the estimated phononic background of k c (cf. 
text). Inset: Estimations for the uncertainties of the phonon 
background in the representative example k c at x = 0.125 
(dotted lines). 



perature, however, is only observed along the c-axis and 
is the signature of strictly one-dimensional magnetic heat 
conduction along the ladders in the materiali^ In the fol- 
lowing we concentrate on the evolution of n c at finite Zn 
doping levels x > which is depicted in Fig. Q]b). Upon 
Zn-doping, n c is progressively suppressed and exhibits 
less pronounced peaks. The magnetic peak simultane- 
ously broadens and apparently shifts towards higher tem- 
perature since the peak structure evolves into a mono- 
tonic increase at high temperature. 

Qualitatively, the doping induced suppression of both 
phononic and magnetic transport channels of k c is a nat- 
ural consequence of the induced structural and magnetic 
defects: different constituents on the same lattice site (Cu 



and Zn) must lead to enhanced scattering of phonons. 
Similarly, since the Zn 2+ ions create nonmagnetic sites 
in a S = 1/2 spin ladder built of predominantly Cu 2+ 
ions, magnetic excitations scatter on these sites. The up- 
per limits for Z mag thus created correspond to the mean 
distances between a pair of Zn-ions along a ladder, i.e. 
rfzn-Zn ~ 374 A, 187 A, 94 A and 62 A. Apparently, even 
the small Zn-contents discussed here lead to an efficient 
suppression of K mag with respect to the undoped mate- 
rial. We can therefore conclude without further analy- 
sis that / mag of pure Sri4Cu2404i must be significantly 
larger than the highest cfen-Zn ~ 374 A. 

Our aim is now to analyze the magnetic heat conduc- 
tivity K mag of the ladders in more detail; the phononic 
contribution K p h to k c must therefore be estimated and 
subtracted from k c . The estimation of K p h is based 
on the fact that due to the large spin gap of order 
~ 400 K^2±. K mag is negligibly small at T < 50 K and 
hence K p h ~ k c in this temperature range. Therefore, the 
usual procedur e) 1 ! 2 ! 5 is to fit k c at T < 50 K with a Debye 
model— and then extrapolate the fit function towards 
higher temperature. In our analysis of K p h in the Zn- 
doped compounds we impose the additional constraint 
of identical fit parameters for all doping levels provided 
that they are not related to phonon-defect scattering. 
The coefficient for phonon-defect scattering is allowed to 
vary with increasing Zn-contentS and the results of this 
procedure are shown in Fig. \T\b) as solid lines. Such an 
extrapolation of K p h is certainly not exempt from uncer- 
tainty. We take these into account by estimating positive 
and negative deviations from Kp^^ A representative ex- 
ample of such deviation is shown in the inset of Fig. ^ for 
x = 0.125. 

Fig-Et displays the resulting K mag of the four Zn-doped 
compounds as a function of temperature. The peak-like 
shape apparent in the curve at x = 0.125 reflects the typ- 
ical temperature dependence of the thermal conductivity 
of a system with a moderate number of impurities: at low 
temperature, where scattering rates are only weakly tem- 
perature dependent, the strong increase in the number of 
quasiparticles generates an increasing k. The decrease of 
k at intermediate temperature (180 K < T < 250 K) and 
the almost constant to slightly increasing k at further ris- 
ing temperature signal a strong increase of the scattering 
rate and a weakening or even saturation of this increase, 
respectively. The relevant scattering processes for K mag 
in the undoped compound Sri4Cu2404i are scattering on 
defects at low temperature and magnon-hole scattering 
at intermediate and high temperatures m& In the case of 
our weakly doped samples with a Zn-concentration < 3% 
it is reasonable to assume that this situation is in prin- 
ciple preserved with the particularity, however, that the 
defects relevant for the magnon scattering are mainly de- 
fined by the Zn-dopants. In this scenario one expects 
the temperature independent magnon-defect scattering 
to increasingly dominate over the temperature depen- 
dent magnon-hole scattering, when the Zn-doping, i.e. 
the density of defects is further increased. K mag evolves 
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FIG. 2: K msg (a) and i mag (b) of Sri4Cu24- x Zn x 04i at 
x = 0.125, 0.25, 0.5, 0.75 and x = (Inset). In (b) data 
points below 75 K are omitted, because the experimental er- 
rors are to large in this temperature range. The grey shaded 
areas in (a) and (b) display the experimental uncertainty of 
Atmag and Z mag resulting from the uncertainties of jc p h for the 
representative case x = 0.125. The solid lines in (c) represent 
the estimated i^ ag (cf. text). The dotted lines mark the mean 
distance of Zn-ions in the ladders. 



accordingly upon doping, since K mag decreases and the 
peak-like structure vanishes at x > 0.25. 

We proceed now with a quantitative analysis of K mag , 
where our main interest lies in a comparison of the mean 
distance of Zn-ions cizn-Zn and £ mag extracted from the 
data. The analysis is based on a kinetic description of the 



heat transport in the ladders presented earlier in Ref. 0, 

_ 3A/7 mag /" £max exp(e/fc B r) 2 
Nnag " nKk B T* J A (exp( e /fc B T) + 3)= 6 ^ W 

Here, N is the number of ladders per unit area; l mag , 
A and e max « 200 meV— are the mean free path, the 
gap and the band maximum of the spin excitations re- 
spectively. Presuming that the spin gap£ A = 396 K 
of undoped Sri4Cu2404i does not significantly change at 
these light Zn-doping levels^! Z mag can directly be calcu- 
lated from our K mag data as shown in Fig. Et>- 

Similar to Z mag (T) at x = (cf. inset of Fig. El, ? mag 
of the Zn-doped samples clearly decreases upon heat- 
ing. The decrease itself, however, is much weaker for 
the doped samples, and apparently weakens further with 
increasing Zn-content. In Ref. H it has been argued that 
in the undoped compound the temperature dependence 
of i roa g is intimately related to temperature dependent 
magnon-hole scattering which arises due to changes of 
hole mobility upon charge ordering: at T < 100 K, i.e. 
deep in the charge ordered state, magnon hole scatter- 
ing is unimportant. It becomes increasingly dominant 
at higher temperature and reaches its full strength at 
T > 240 K, i.e. above the charge ordering temperature. 
Z mag is thereby reduced from l Q w 3000 A at T < 100 K 
to around 60 A at 300 K. 

Following Matthiesen's rule it is possible to separate 
tag = ^cT 1 + contributions due to scattering off 

static defects, lo, and off mobile holes, The observed 
reduction of Z mag upon Zn-doping is indeed consistent 
with a gradually reduced Iq due to scattering of magnons 
on Zn-ions and only a slight change in 1^. In order to 
illustrate this, we estimate the mean free path of the 
Zn-doped samples under simplified conditions, i.e. we 
assume that firstly lo = (izn-Zn and secondly that lh{T) 
is the same as in the non-doped material as published 
in Ref. 0| hence Z*~ ag = dzn-Zn -1 + l^ 1 .^ The results 
for Z^ ag are displayed in Fig. Et- For T > 100 K there 
is a striking similarity in both magnitude and tempera- 
ture dependence between the results of the experimental 
Z mag and our estimated Z^ ag values. We should mention 
the unavoidable uncertainties in ^ mag which are large at 
T < 100 K, since the ratio K mag /Kph is small here and 
errors in the estimation of n p h have large effects on the 
estimations of K mag and i mag (cf. Fig.EK andEb)- A reli- 
able determination of the experimental / mag can therefore 
only be expected at higher temperature. 

However, at 100 K < T < 125 K such errors be- 
come reasonably small and simultaneously Z mag ~ dzn-Zn 
should hold. (Note that in this temperature range the 
deviation of Z^ ag from Iq = dzn-Zn is smaller than 10%, 
as can be inferred from the dotted lines in Fig.Et-) We 
therefore test our experimental Z mag from Fig.Eb for such 
a linear scaling at a representative T — 120 K in Fig. 01 
As is evident from the figure, i mag roughly scales with 
^Zn-Zn and a linear fit to the data points yields a slope 
of 1.17 ± 0.25, which is close to unity. This confirms un- 
ambiguously that the values for i mag obtained from the 
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FIG. 3: i mag as a function of dzn-Zn- Solid line: linear fit line 
through the origin. The error bars arise due to uncertainties 
in determining re p h and have been obtained by using different 
estimates for K p h (cf. inset of Fig. and Fig. 



kinetic model are in fair agreement with real scattering 
lengths in the material. 



IV. DISCUSSION 

The consequences of the analysis outlined above are 
evident. It is natural to conclude that the kinetic model 
yields realistic numbers for Z mag , not only in the Zn-doped 
material discussed above but also in Sri4Cu2404i and 
La5CagCu2404i where it yields the aforementioned sur- 
prising high values of Z mag ~ 3000 A at low temperature. 
Therefore, the rather small values for / mag , which were 
obtained under the assumption of a finite thermal Drude 
weight 1 ^ become less realistic. This is further reinforced 
by recent theoretical investigations which suggest a van- 
ishing thermal Drude weigh1*i2*2ii. 

We point out that on a microscopic level the mech- 



anism of magnon-scattering off nonmagnetic Zn 2+ ions 
and its relevance for magnetic heat transport remains 
unaddressed. Qualitatively this scattering appears to 
be very similar to the 2D case of La2Cu04, where 
( mag w dzn-Zn has been observed as well (with dzn-Zn 
being the mean unidirectional distance of Zn-impurities 
in the CuOVplanes)^ In this regard it seems to be highly 
intriguing to investigate the heat transport of Zn-doped 
spin chain materials such as SrCuC>2 and S^CuOa, since 
Zn-impurities must lead to complete chain interruptions 
in these cases which could - compared to two-leg 
ladders and 2D antiferromagnetic planes - result in a 
qualitatively different effect on K mag - 

V. SUMMARY 

Our experimental study of >t mag of Sri4Cu24- x Zn x 04i 
shows that Zn-doping in this spin ladder compound 
strongly suppresses K mag with a progressive decrease 
upon increasing Zn-content. Since the Zn-ions create de- 
fects with a well defined distance the data provide a clear- 
cut and model-independent relation between K mag and 
the mean free path i mag < 374 A. This strongly corrob- 
orates previous results for l mag of pure Sri4Cu24 04i and 
La5CagCu2404i, in which surprisingly large values for 
/ mag have been found on basis of kinetic model-. We fur- 
ther find that the Z mag determined from the present data 
with this kinetic model are roughly equal to the mean dis- 
tance of Zn-ions along the ladders. Magnetic scattering 
lengths in the spin ladder material (Sr, Ca, La)i4Cu2404i 
can therefore conveniently be deduced from experimental 
data for K m ^. 
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For the fit we used 3 ^ Kph = ^ — j Jq j^pr 

with the scattering rate t^ 1 = Auj 4 + v/L + 
BTto 3 exp(-Q D /bT), where u, Q D and v are the phonon 
frequency, Debye temperature and mean phonon velocity. 
The three terms correspond to scattering due to defects, 
crystal boundaries and umklapp processes respectively. A, 
L, B and b are free fit parameters with the constraint that 
L, B and b are identical for all Zn-doped samples. 
The different phononic backgrounds have been obtained 
from fitting data points at T < (53 ± 12) K. 
Despite other Zn-doped low-dimensional quantum spins 
systems 3 — our data on Sri4Cu24- x Zn x 04i do not yield any 
evidence for magnetic order at low temperature. 
Resistivity data of Sri4Cu24- x Zn x 04i indicate a weaken- 
ing of the charge-ordered state upon Zn-doping. An addi- 
tional suppression of Z mag could therefore arise due to an 
increased mobility of holes in the ladders and a thereby 
reduced In- However, this effect appears to be only of mi- 
nor importance, since a clear signature of charge ordering is 
still observed at the lowest doping level x — 0.125, whereas 
Zmag is already strongly suppressed. In any case, the main 
doping- and temperature dependent features of Z mag are 
quite robust against such doping dependencies as long as 
lh S> Iq at low temperature. 
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